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Abstract 

Cluster algebras are a recent topic of study and have been shown to be a 
useful tool to characterize structures in several knowledge helds. An impor¬ 
tant problem is to establish whether or not a given cluster algebra is of hnite 
type. Using the standard dehnition, the problem is infeasible since it uses 
mutations that can lead to an inhnite process. Barot, Geiss and Zelevinsky 
(2006) presented an easier way to verify if a given algebra is of hnite type, by 
testing that all chordless cycles of the graph related to the algebra are cycli¬ 
cally oriented and that there exists a positive quasi-Cartan companion of the 
skew-symmetrizable matrix related to the algebra. We develop an algorithm 
that verihes these conditions and decides whether or not a cluster algebra is 
of hnite type in polynomial time. The second part of the algorithm is used 
to prove that the more general problem to decide if a matrix has a positive 
quasi-Cartan companion is in MV. 
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1. Introduction 

Matrices can be used to represent various structures arising in combina¬ 
torics, including graphs and algebras, such as cluster algebras. The latter 
can be dehned using a directed graph G{B), called quiver, and consequently 
by an adjcency matrix, where rows and columns represent the vertices and 
the positive values at positions {i,j) represent the quantity of edges between 
associated vertices of the graph. 
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In 2002, Sergey Fomin and Andrei Zelevinsky introdnced a commn- 
tative algebra class called cluster algebras. These algebras have a strong 
combinatorial structure. They are tools to study questions of dual canonical 
bases and positivity of semisimple Lie groups. Cluster algebras are defined 
recursively via commutative algebras with a distinct set of generating vari¬ 
ables {cluster variables) grouped into overlapping subsets {clusters) of hxed 
cardinality. 

A basic feature of any cluster algebra class is that both the generators 
and the relationships between them are not given from the start, but are 
produced by an elementary iterative process of seed mutation. This process 
is somewhat counterintuitive, in order to encode a universal phenomenon in 
some way. This may explain the accelerated development of cluster algebra 
theme in areas such as combinatorics, physics, mathematics (especially ge¬ 
ometry), among others as discussed in (9|, llO, ll2| . This algebra can be defined 
using a skew-symmetrizable matrix, as we will see in Section [3l 

The notion of quasi-Cartan matrices was introduced by Barot, Geiss and 
Zelevinsky [l|. A quasi-Cartan matrix is a symmetrizable matrix with all 
entries of the main diagonal equal to 2. The authors show some properties 
of the matrices from the mathematical point of view. A quasi-Cartan com¬ 
panion is a quasi-Cartan matrix associated with skew-symmetrizable matrix 
as we will see in Section [2l 

One can decide whether a cluster algebra is of finite type (has a finite num¬ 
ber of cluster variables) by deciding whether or not the skew-symmetrizable 
matrix is associated with a cyclically oriented graph and has a positive quasi- 


Cartan companion. By the Sylvester criterion [1^, a symmetric matrix is 


positive if all its leading principal submatrices have positive determinant. 

In this paper, we present three algorithms of polynomial time complexity. 
The hrst one can be used to decide whether or not a quasi-Cartan companion 
matrix is positive. This is used as a certificate to prove that the problem of 
setting whether a positive quasi-Cartan companion exists or not belongs to 
the AfV class of problems. The second algorithm can be used to establish 
whether or not there exists a positive quasi-Cartan companion of a skew- 
symmetrizable matrix associated with a cyclically oriented graph. To verify 
that an oriented graph is cyclically oriented, we will use the algorithm and 
results presented in j2|. The last algorithm can be used to decide whether or 
not a cluster algebra is of hnite type. This is used to prove that the problem 
belongs to the V class of problems. For more information about V, AfV and 
AAP-complete classes, see BQ. 


2 






The remainder of the paper is organized as follows: Section [2] presents 
the preliminary concepts; Section [3] dehnes what is a cluster algebra of £- 
nite type and presents some of its properties; In Section 01 we present some 
properties of cyclically oriented graphs and the general idea of the polyno¬ 
mial algorithm for its determination presented is in ; Section |5] presents 
a polynomial algorithm for determining if a skew-symmetrizable matrix has 
a positive quasi-Cartan companion when the associated graph is cyclically 
oriented; Section [6] shows that deciding whether a cluster algebra is of hnite 
type is a problem belonging to the class of polynomial problems. Finally, 
Section [7] closes the paper with concluding remarks and future work. 

2. Preliminaries 

Let n be a positive integer, A,B,C G M„(Z) and D G M„(Q). A matrix 
A is symmetric if A = A^, where A^ is the transpose of A. A matrix C is 
symmetrizable ii DxC is symmetric for some diagonal matrix D with positive 
diagonal entries. In this case, the matrix D x C is the symmetrization or 
symmetrized version of C and the matrix D is the symmetrizer of C. Note 
that this dehnition is equivalent to the one given in |3| and the matrix D is 
taken over Q for the sake of simplicity. A matrix C is symmetric by signs if 
for all i, j G {1,..., n}, with i ^ j, we have Cij = Cji = 0 or Cij ■ Cji > 0. 

All symmetric matrix are symmetrizable and that all symmetrizable ma¬ 
trices are symmetric by signs. 

A matrix A is skew-symmetric if A^ = —A. Observe that the values 
of the main diagonal are null. A matrix B is skew-symmetrizable if there 
exists a diagonal matrix D with positive entries such that D x B is a. skew- 
symmetric matrix. In this case, the matrix DxB is the skew-symmetrization 
or skew-symmetrized of B and the matrix D is the skew-symmetrizer of B. 
A matrix B is skew-symmetric by signs if for all i,j G {1,... ,n} we have 
bii = 0 and if i ^ j, then bij = bji = 0 or hij ■ hji < 0. 

Also observe that all skew-symmetric matrix are skew-symmetrizable and 
that all skew-symmetrizable matrices are skew-symmetric by signs. 

A generalized Cartan matrix is a symmetric matrix, where its main diag¬ 
onal entries are equal to 2 and its other values are non-positives. A Cartan 
matrix is a symmetrizable matrix with a positive dehnite symmetrized ma¬ 
trix. Cartan matrices were hrst introduced by the French mathematician 
Elie Cartan. In fact, Cartan matrices, in the context of Lie algebras, were 
hrst investigated by Wilhelm Killing, whereas the Killing form is due to 
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Cartan. Positive Cartan matrices represent the basis of the Cartan-Killing 
classihcation, see (sf. 

A symmetrizable matrix is quasi-Cartan if all the entries on its main 
diagonal equal to 2. For a skew-symmetrizable matrix B, we will refer to 
a quasi-Cartan matrix C with \cij\ = |%| for all i 7 ^ j as a quasi-Cartan 
companion of B. 

Given a skew-symmetrizable matrix B, we associate G{B) to an oriented 
graph with vertices { 1 , 2 ,, n} and edges {i,j) for each bij > 0 , with i,j e 




Let G be a simple graph, with edge set E and vertex set V. Let n 
be the number of vertices in V and m the number of edges in E. Cyclical 
orientable graphs were introduced in 2006 by Barot, Geiss and Zelevinsky [l|. 


A simple path is a hnite sequence of vertices {vi,V 2 , ■ ■ ■ ,Vt) such that 
{vi, Vj+i) G E for i G — 1} and no vertex is repeated in the sequence, 

that is, Vi 7 ^ Vj, for i,j G {1,... ,f} and i ^ j- A cycle is a simple path 
(^ 1 ,^ 2 , ■ ■ ■ ,Vt) such that {vt,vi) G E. We denote a cycle with t vertices by 
Gto A chord of a cycle is an edge between two vertices of the cycle, that is 
not part of it. A cycle without chord is a chordless cycle. 

An orientation of a graph G{B) consists of assigning an order to the end¬ 
points of each of its edges. The orientation of a cycle {vi,V 2 ■ ■ ■, Vt) is cyclic 
if it receives the orientations (ui, ^ 2 ),..., (ut-i, ft), (ft, Ui) or the opposite. 

An oriented graph G{B) is cyclically oriented if any chordless cycle in 
G{B) is cyclic. It is cyclically orientable if it admits an orientation in which 
G{B) is cyclically oriented. 

3. Cluster algebras of finite type 

Gluster algebras are associated with the mutation-equivalence classes of 
skew-symmetrizable matrices. Recall from that, for each matrix index k, 
the mutation in direction k transforms a skew-symmetrizable matrix B into 
another skew-symmetrizable matrix B’ = fik{B), whose entries are given by 



bij, 

bij + sgn{bik) ■ [kk ■ bkj]+, 


a i = k or j = k] 
otherwise. 


( 1 ) 


Observe that our definition of a cycle, as in Q, does not repeat the first vertex at the 
end of the sequence as usually done by other authors. 
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where we use the notation [a:]_|_ = max(a:, 0 ), with the convention sgn{0) = 0 
and sgn{x) = a:/|a:|. One can easily check that /i^ is involutive, implying that 
the repeated mutations in all directions give rise to the mutation-equivalence 
relation on skew-symmetrizable matrices. Each of these matrices dehne an 
exchange rule. 

For an integer n, a cluster algebra of dimension n is a commutative ring 
(unitary) without zero divisors, generated in the center of a hxed held T, 
for a set (possibly inhnite) of cluster variables. These variables are not hxed 
arbitrarily. The set of cluster variables is the union (not disjoint) of subsets of 
n elements called clusters, which are related by the mutation change. For any 
cluster X and any cluster variable a: G <T, there is another cluster obtained by 
substitution of the variable x by another x' related by a binomial relationship 
as follows: x-x' = M 1 + M 2 , where Mi and M 2 are dehned using an exchange 
rule matrix. All cluster variables are recursively obtained in this way from 
an initial seed (a cluster with an exchange rule matrix). This process is a 
seed mutation. For more information, see j^. 

Cluster algebras of hnite type are those that have a hnite number of 
cluster variables. They coincide with the famous Cartan-Killing classihcation 
of semisimple Lie algebras. 

The following criterion for deciding whether a skew-symmetrizable matrix 
corresponds to a cluster algebra of finite type is presented in [l|. 

Theorem 1 (Barot, Geiss and Zelevinsky |[l|). Given a cluster algebra 
A{B) associated with a skew-symmetrizable matrix B. LetS be the mutation- 
equivalence class of B. The following are equivalent: 


(a) The cluster algebra A{B) is of finite type. 

(b) S contains a matrix B' such that the Cartan matrix C with off-diagonal 
entries Cij = —\hT\ is positive. 

(c) For every B' E S and all i 7^ j, we have \bT ■ 6'J < 3. 

(d) Every chordless cycle in G{B) is cyclically oriented, and B has a pos¬ 
itive quasi-Gartan companion. 


Furthermore, the Gartan-Killing type of the Gartan matrix C in (b) is uniquely 
determined by S. 
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Mutations are hard to control, so each of the conditions (b) and 


in 


Theorem [T] is hard to check in general, since there may be a very large (pos¬ 


sibly inhnite) number of matrices in S. On the other hand, the condition (d) 
leads us to a polynomial algorithm as we see in Section |6l As with condi¬ 
tion (d), we divide the verihcation into two parts: Section 0] presents the 


ideas of the algorithm to decide whether a graph G{B) is cyclically oriented 
(that was proposed by us in j^) and in Section Owe present two polynomial 
algorithms that, together, decide whether B has a positive quasi-Cartan com¬ 
panion when G{B) is cyclically oriented. 


4. Cyclically oriented 


In this section, we will present some properties of chordless cycles, cycli¬ 
cally orientable and cyclically oriented graphs. 

Some properties of cyclically orientable graphs are given by Barot, Geiss 
and Zelevinsky [l|. New characterizations were obtained by Speyer [l^ and 
Gurvich 0 in later works, which enabled the development of algorithms for 
the recognition of cyclically orientable graphs. 

The proposition below is used in preprocessing of our algorithm to deter¬ 
mine whether a graph is cyclically oriented. 


Proposition 1 (Speyer jl4l|i. If G is a cyclically orientable graph with n 
vertices, then G has at most 2 ■ n — 3 edges. 


Now we will discuss two-connected components. They are important 
since any chordless cycle is contained in exactly one of the components. To 
calculate them, we can use an algorithm based on the ideas of Tarjan 0 
and Szwarchter [l^, that has time complexity which is an application 

of depth-first search (DFS), described in Cormen et ah j^. 

A graph G is connected when there exists a path between each pair of 
vertices of G, otherwise G is disconnected. A connected component of a 
graph G is a maximal connected subgraph of G. A graph is two-connected 
if it is connected and the elimination of at least two vertices is necessary to 
disconnect it. A single edge is a two-connected graph. 

The following theorem deals with cyclic orientability in two-connected 
components. The idea is used in our algorithms. 


Theorem 2 (Speyer jl4l|f . A graph G is cyclically orientable if and only 
if all of its two-connected components also are. A two-connected graph is 
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cyclically orientable if and only if it is either a cycle, a single edge, or of 
the form G' U Ct, where G' is a cyclically orientable graph, Gt is a cycle and 
G' and Gt meet along a single edge. Moreover, if G = G' U Gt is any such 
decomposition of G into a cycle and a subgraph meeting along a single edge, 
then G is cyclically orientable if and only if G' is. 

4 . 1 . Algorithm to determine whether a graph is cyclically oriented 

Based on the results given by Speyer [l^, Dias and Castonguay in 
proposed a polynomial algorithm to verify if a graph is cyclically oriented. 
The algorithm can be slightly modihed to verify whether or not a given 
oriented graph is cyclically oriented. Moreover, in the positive case, they 
show that the algorithm returns all chordless cycles of the graph. 

We now highlight the ideas of the algorithm. It is based on the analysis of 
each two-connected component of a given graph. Following the idea presented 
by Speyer jl^, the algorithm identihes chordless cycles induced by cycle Gt 
given by Theorem |2l This is done in order to reduce each two-connected 
component to a unique cycle. 

Initially, the algorithm verihes if the given graph satishes Proposition [H 
that is, if the graph has 2 • n — 3 edges. If not, it returns the decision 
NO. Next, it hnds all two-connected components and then verihes if each 
component satishes Proposition [H If not satished, it returns NO. 

After that, the algorithm stores in a queue Q all vertices of degree two 
for each two-connected component. Observe that all vertices will be added 
to Q at most once and if G is cyclically orientable then all vertices will be 
added. 

The algorithm tries, starting with the vertices of queue Q, to hnd and 
eliminate paths (cycles) in order to reduce the initial two-connected compo¬ 
nent to a cycle. If it is able to do so, the component is cyclically oriented. If 
not, it returns NO. This will continue for all components. In the hnal process, 
the given graph will be classihed as cyclically orientable if all two-connected 
components receive a cyclically orientable classihcation; otherwise, the graph 
is classihed as not cyclically orientable. Therefore, it determines in 0{n^) 
time whether this graph is cyclically orientable and, if so, it returns the set 
of all chordless cycles S of G. Here, we will consider the set S' as a stack. 
Observe that there is at most n chordless cycles in the whole graph. 

Dias and Castonguay provided an algorithm to verify whether a non- 
oriented graph is cyclically orientable. Observe that this paper consider 
soriented graphs, therefore it is necessary to make small modihcations in 
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order that the algorithm can be used in Algorithm [31 Consider the modihed 
algorithm called ChordlessCyclesCOd{G). If a graph is cyclically oriented 
it returns the set S of all chordless cycles. For more information about the 
algorithm, see j^. 


5. Positive quasi-Cartan companion 

As we saw in Theorem [H Barot, Geiss and Zelevinsky [l| show that for a 
cluster algebra of hnite type we have that the associated skew-symmetrizable 
matrix B must have a positi ve q uasi-Cartan companion C. 

By the Sylvester criterion 1^, a symmetric matrix is positive if all leading 
principal submatrices have positive determinant. A leading principal subma¬ 
trix is obtained by iteratively removing the last row and the last column 
of the matrix. Using j^, a symmetrizable matrix is positive if all leading 
principal submatrices have positive determinant. 

Dias, Castonguay and Dourado presented Algorithm [1] with time com¬ 
plexity 6{n^) to decide whether the given matrix C is positive. This algo¬ 
rithm is used as a veriher for the general problem of deciding if there exists 
a positive quasi-Cartan companion. Thus, the problem belongs to the MV 
class. 


Algorithm 1: IsPositive{C) 

Input: A n X n symmetrizable matrix C. 

Output: The response is if the matrix is positive or not. 

1 foreach leading principal submatrix C of C do 

2 if (det{C') <0) then 

3 |_ return NO 

4 return YES 


Verifying all quasi-Cartan companions is exponential since if i? is a n x n 
skew-symmetrizable matrix there are 2™ matrices C which are quasi-Cartan 
companions of B. Note that a quasi-Cartan companion of B is specified by 
choosing the signs of its off-diagonal matrix entries, with the only requirement 
being that sgrifcij) = sgn{cji) and \cij\ = |6jj|, for i ^ j. Dias, Castonguay 
and Dourado [6| showed that if C is symmetric according to the signs, then 
C is symmetrizable with the same symmetrizer of B. 








Thus, it is not efficient to test all possible signs for a positive quasi-Cartan 
companion. However, Barot, Geiss and Zelevinsky [l| showed a sign condition 
which simplihes the test of the existence of a quasi-Cartan companion without 
having to assign all possible signs. 

Definition 1 (Sign condition — Barot, Geiss and Zelevinsky |[l|). Let 

B be skew-symmetrizable matrix. A quasi-Cartan companion C of B satis¬ 
fies the sign condition if, for every chordless cycle Ct in G{B), the product 
n (“Cjj) over all edges of Ct is negative. 

(ij)GCt 


In our case, the skew-symmetrizable matrix always has at least one quasi- 
Cartan companion that satishes the sign condition. 

Proposition 2 (Barot, Geiss and Zelevinsky |[l|). Let B be a skew-symmetrizable 
matrix. If G{B) is cyclically oriented, then B has a quasi-Cartan companion 
(not necessarily positive) satisfying the sign condition; furthermore, such a 
quasi-Cartan companion is unique up to simultaneous sign changes in rows 
and columns. 

The next proposition shows that it is enough to check the quasi-Cartan 
companion satisfying the sign condition. 

Proposition 3 (Barot, Geiss and Zelevinsky jit]). To be positive, a quasi- 
Cartan companion C of a skew-symmetrizable matrix B must satisfy the sign 
condition. 

Observe that a quasi-Cartan matrix that satishes the sign condition is 
not necessarily positive. 

The following proposition shows that, for our purpose, it is enough to 
verify the positivity of one quasi-Cartan companion that satishes the sign 
condition. The proposed algorithm to verify if a skew-symmetrizable matrix 
has a positive quasi-Cartan matrix is based on hnding a quasi-Cartan matrix 
that satishes the sign condition. 

Proposition 4. Let C be a quasi-Cartan companion of a skew-symmetrizable 
matrix B which satisfies the sign condition. If G{B) is cyclically oriented, 
then B has a positive quasi-Cartan companion if and only if C is positive. 


9 


Proof. Suppose that G{B) is cyclically oriented. Clearly, if C is positive, 
then B has a positive quasi-Cartan companion. On the other hand, if B 
has a positive quasi-Cartan companion, say C, then by Proposition [3l C 
satishes the sign condition. By Proposition [2l C is obtained from C by 
simultaneous sign changes in rows and columns. Therefore, there exists a 
diagonal matrix X = {xij) with xu G {—1,1} such that C = XC'X. Since 
det{C) = det{X) ■ det{C') ■ det{X) and det{X) G {—1,1}, we have that C is 
positive. □ 

Next we present an algorithm that dehnes an attribution of sign for a 
quasi-Cartan companion of B. In the algorithm, the obtained matrix C 
will satisfy the sign condition. Later, we verify the positivity of C. By 
Proposition 01 the matrix B has a positive quasi-Cartan companion if and 
only if C is positive. 

Observe that the stack S contains all chordless cycles obtained by the Al¬ 
gorithm ChordlessCyclesCOd{G) and the set T contains all two-connected 
components of G that are single edges. 

Based on , one can see that the cardinality of S is at most n. Therefore, 
the time complexity of Algorithm [2] is 0{n‘^), due Line [191 The rest of the 
algorithm is 0{n‘^). To prove the correctness of Algorithm [ 2 l we hrst show 
that the function sgn{) is well dehned. 

Proposition 5. The function sign is well defined, that is, sgn{edge) G {1, —1} 
for all edge G E. 

Proof. The line numbers in this proof refer to lines of Algorithm [2l First, we 
prove that in Line [151 L 7 ^ 0. We can assume that G is a two-connected graph, 
since the algorithm acts on each two-connected component. The hrst cycle 
is a cycle without a dehned edge. Therefore, the hrst edge does not satisfy 
the condition in Line [91 and enters on Line 021 (/ = 1). For any other cycles 
Gt, we have G = G' U Gt such that only edges of G' are dehned. Therefore, 
the hrst or second edge of Gt is not dehned (J = 1 or J = 2 ). 

Since the graph G is two-connected, for each connected component we 
have a single edge or the edge is part of a cycle. If G is a single edge, then 
it receives 1 at Line [51 If not, the edge receives 1 at Line 091 or “-prod” at 
Line 001 Therefore, sign is well-dehned. □ 
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Algorithm 2: PositiveCompanion{G, B, S) 


Input; A cyclically oriented graph G, a skew-symmetrizable matrix B 

associated to G, a stack S of chordless cycles and a set T of single 
edges (two-connected components of G). 

Output: The response whether or not there exists a positive quasi-Cartan 
companion. 

1 initialize sgn{xi, Xj) ■(— 0 for all i,j G {1, 2,... , n} 

2 initialize sgn{edge) G- 1 for all single edges in T 

3 while {S / 0) do 

4 remove an element c of S' /* c = (xi,..., Xt) */ 

5 Xt+l G- Xi 

6 I 0 

7 prod G- 1 

8 foreach i G t} do 

9 if {sgn{xi,Xi+i) 7 ^ 0) then 

10 I prod G- prod ■ sgn{xi, Xj+i) 

11 else 

12 if (/ = 0 ) then 

13 I I i 

14 else 

15 |_ sgn{xi,Xi+i) G- 1 

16 sgn{xi,xi+i) < - prod 

17 initialize Cij G- \bij \ ■ sgn(xi,Xj) for all i,j G {1,2,... ,n} 

18 initialize cu G- 2 for all i G {1, 2,... , n} 

19 if {IsPositive{G)) then 

20 I return YES 

21 else 

22 1^ return NO 


We show next that the construction of C gives a quasi-Cartan companion 
that satisfies the sign condition. 

Theorem 3. The n x n matrix C defined below satisfies the sign condition. 



2 , 

sgn{xi,Xj) ■ \bij\, 


ifi =j 

otherwise. 
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Proof. Recall that S is composed of all chordless cycles of G. It follows from 
Proposition O that the product of each cycle is equal to —prodf. □ 


Theorem 4. The Algorithmic is correct. 

Proof. Follows from j^, Theorem |3] and Proposition |H □ 


6 . Recognizing cluster algebras of finite type is in T* 

Since the algoritm ChordlessCyclesCOd{G) checks whether or not a 
graph is cyclically oriented and the algorithm PositiveGompanion{G, B, S) 
verihes whether or not the matrix has a positive quasi-Cartan companion, 
both in polynomial time complexity, it follows from Theorem [T] that decid¬ 
ing whether or not a cluster algebra is of hnite type belongs to the class of 
polynomial problems {V). 


Algorithm 3: GlusterAlgebraFiniteType{B,G) 

Input: A skew-symmetrizable matrix B that defines a cluster algebra 
A{B), and the respective graph G. 

Output: The response that the cluster algebra A{B) is of finite type or 
not. 

1 COd, S'^ ChordlessCyclesCOd(G') 

2 if (COd = YES) then 

3 I PositiveCompanionC GjRS) 

4 else 

5 return NO 


Based on the above algorithms, we have that the time complexity of 
Algorithm [3] is 0{n^). The correctness of Algorithm [3] follows directly from 
the correctness of Algorithms [1] and [2l 

7. Conclusions and future work 

In this paper, we proposed two polynomial algorithms that together prove 
that deciding if a cluster algebra is of finite type belongs to the class of poly¬ 
nomial problems {V). The first checks whether or not the associated graph is 
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cyclically orientable and, if so, the second verifies whether there exists a posi¬ 
tive qnasi-Cartan companion of B, which represents a clnster algebra. These 
two criteria can be evalnated for the determination of hniteness of a clnster 
algebra due to Theorem [1] presented by Barot, Geiss and Zelevinsky jl|. 

The result shown is important because the use of an efficient algorithm to 
determine if a cluster algebra is of hnite type will facilitate research in many 
areas of the application of cluster algebras. Another interesting solution of 
the problem was obtained by A. Seven in [l^, in which it is verihed whether 
a cluster algebra is of hnite type in terms of “forbidden minors” of B. 

Importantly, we conjecture that, in general, the problem of hnding a 
positive qnasi-Cartan companion for a skew-symmetrizable matrix B is an 
AAP-complete problem. This problem will be considered in future work. 
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